1 Call Libraries and Set Options

library(tidycensus) # get tidy data from the Census

library(tidyverse) # tidyverse including ggplot

options(tigris_use_cache = TRUE) # options

library(mapview) # interactive maps

2 How Do I Know What Geographies or Variables I Want?

https://walkerke.github.io/tidycensus/articles/basic-usage.html

library(tidycensus)

v17 <- load_variables(2017, "acs5", cache = TRUE)

# View(v17)

library(DT)

datatable(v17,
          filter = 'top',
          options = list(dom = 'Bfrtip', 
                         buttons = c('copy', 
                                     'csv', 
                                     'excel', 
                                     'pdf', 
                                     'print')))

3 Census API Key

census_api_key("...") # Census API key

4 Get Data

Washtenaw <- get_acs(state = "MI", 
                     county = "Washtenaw", 
                     geography = "tract", 
                     variables = "B19013_001", 
                     geometry = TRUE)

5 Look at Data

head(Washtenaw) # look at data

## Simple feature collection with 6 features and 5 fields
## geometry type:  MULTIPOLYGON
## dimension:      XY
## bbox:           xmin: -83.75863 ymin: 42.26068 xmax: -83.72033 ymax: 42.29117
## epsg (SRID):    4269
## proj4string:    +proj=longlat +datum=NAD83 +no_defs
##         GEOID                                          NAME   variable
## 1 26161400100 Census Tract 4001, Washtenaw County, Michigan B19013_001
## 2 26161400200 Census Tract 4002, Washtenaw County, Michigan B19013_001
## 3 26161400300 Census Tract 4003, Washtenaw County, Michigan B19013_001
## 4 26161400400 Census Tract 4004, Washtenaw County, Michigan B19013_001
## 5 26161400500 Census Tract 4005, Washtenaw County, Michigan B19013_001
## 6 26161400600 Census Tract 4006, Washtenaw County, Michigan B19013_001
##   estimate   moe                       geometry
## 1    34625 10975 MULTIPOLYGON (((-83.75217 4...
## 2    38403 16912 MULTIPOLYGON (((-83.7431 42...
## 3    34280  9044 MULTIPOLYGON (((-83.74315 4...
## 4   109338 15489 MULTIPOLYGON (((-83.74047 4...
## 5    22365  8032 MULTIPOLYGON (((-83.75064 4...
## 6    59551 14879 MULTIPOLYGON (((-83.75863 4...

6 Make a Map

Washtenaw %>% # data
  ggplot(aes(fill=estimate)) + # ggplot
  geom_sf() + # geom for shapefiles
  labs(title = "Quick Map: Washtenaw County",
       subtitle = "Median Family Income") +
  scale_fill_viridis_c() # viridis colors

7 Interactive Version

mapview(Washtenaw, 
        zcol = "estimate", 
        legend = TRUE, 
        layer.name="Median Family Income")

8 Try a Different Theme and Colors

https://github.com/gadenbuie/ggpomological

library(ggpomological)

Washtenaw %>% # data
  ggplot(aes(fill=estimate)) + # ggplot
  geom_sf() + # geom for shapefiles
  labs(title = "Quick Map: Washtenaw County",
       subtitle = "Median Family Income") +
  scale_fill_distiller(palette = "Greens") +
  theme_pomological_fancy()

9 ggthemes

https://jrnold.github.io/ggthemes/

library(ggthemes)

Washtenaw %>% # data
  ggplot(aes(fill=estimate)) + # ggplot
  geom_sf() + # geom for shapefiles
  labs(title = "Quick Map: Washtenaw County",
       subtitle = "Median Family Income") +
  theme_wsj() +
  theme(title = element_text(size = rel(.75)),
        axis.text = element_text(size = rel(.5)),
        legend.text = element_text(size = rel(.5)))

Washtenaw %>% # data
  ggplot(aes(fill=estimate)) + # ggplot
  geom_sf() + # geom for shapefiles
  labs(title = "Quick Map: Washtenaw County",
       subtitle = "Median Family Income") +
  theme_economist() +
  scale_fill_distiller(palette = "greens",
                       direction = 1) +
  theme(title = element_text(size = rel(.75)),
        axis.text = element_text(size = rel(.5)),
        legend.text = element_text(size = rel(.5)))

Washtenaw %>% # data
  ggplot(aes(fill=estimate)) + # ggplot
  geom_sf() + # geom for shapefiles
  labs(title = "Quick Map: Washtenaw County",
       subtitle = "Median Family Income") +
  theme_fivethirtyeight() +
  theme(title = element_text(size = rel(.75)),
        axis.text = element_text(size = rel(.5)),
        legend.text = element_text(size = rel(.5)))

library(ggthemes) 

Washtenaw %>% # data
  ggplot(aes(fill=estimate)) + # ggplot
  geom_sf(alpha=.75) + # geom for shapefiles
  labs(title = "Quick Map: Washtenaw County",
       subtitle = "Median Family Income") +
  scale_fill_gradient(low = "green",
                      high = "darkgreen") + 
  theme_hc(bgcolor = "darkunica") +
  theme(title = element_text(size = rel(.75),
                             color = "green"),
        axis.text = element_text(size = rel(.5),
                                 color = "green"),
        legend.text = element_text(size = rel(.5),
                                   color = "green"))